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Abstract. 

The main objective of the Penn State/Toruh Centre for Astronomy Search for Planets around 
Evolved Stars is the detection of planetary systems around massive, evolved stars. We are also 
interested in the evolution of these systems on stellar evolution timescales. In this paper we 
present our approach to determine the basic physical parameters of our targets GK-giants. We 
also discuss the stellar activity indicators used in our survey: line bisector and curvature, and 
Ha variability. 
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1. Introduction 

Proper interpretation of the results from precision RV studies of GK-giants requires a 
detailed knowledge of their physical parameters. Effective temperatures and gravitational 
accelerations are needed to obtain luminosities, and, with the additional knowledge of 
metallicities, estimates of stellar masses and ages can be derived by means of the isochrone 
fitting. Together with estimates of stellar radii and rotation periods, these data allow one 
to address the influence of stellar surface inhomogeneities (spots) on the observed RV 
variations. 

All alternative sources of RV variations in GK-giants have to be ruled out before 
substellar companion interpretation becomes acceptable. Unfortunately, the long period 
variations, if present, cannot usually be studied using data other than the existing pho- 
tometry (usually of moderate quality), or the data collected from the RV survey. There- 
fore, the detailed activity discussion is usually based on the indicators defined on the 
basis of the same spectra as those used for the RV measurements. In this paper, we illus- 
trate our methodology of the determination of stellar parameters and activity analysis 
based on the case of the K-giant PSU-TCfA 18, a potential planet hosting star. 

2. Observations 

Observations were made with the Hobby-Eberly Telescope (HET) and the High Reso- 
lution Spectrograph (HRS). The HRS was used in the R=60,000 resolution mode with a 
gas cell (I2) inserted into the optical path, and it was fed with a 2 arcsec fiber. Typically, 
the signal-to-noise ratio per resolution element (at 594 nm) was ~200 for the stellar 
spectra taken with the gas cell, and ^250 for the templates. 
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3. Basic stellar parameters 

The atmospheric parameters of the program stars were obtained with the spectroscopic 
method (Takeda et al. 2005a, b) , which is based on analysis of Fe I and Fe II lines and 
relies on conditions resulting from the assumption of the LTE. Typically, over 200 Fel 
and about 25 Fell lines were measured for every star. 

We tested a reliability of our determinations with the Takeda et al. (2005a) TGVIT 
code by applying it to 8 stars, for which the parameters have been published by Butler 
et al. (2006). A comparison of the results shows that T e // values agree to within 49 K, 
and that the same is true for log g to within 0.11 dex and [Fe/H] to within 0.11 dex, 
respectively. 

Stellar masses were derived by comparing the positions of stars in the HR diagram with 
the theoretical evolutionary tracks of Girardi et al. (2000) and Salasnich et al. (2000) for 
a given metallicity. For stars for which the parallax determinations are precise enough, 
the metallicity may introduce a significant uncertainty in mass because of the choice of 
an evolutionary track. We assume that for an average red giant with a known parallax, 
the mass may be estimated to within 0.3M Q . We also note that for stars in the red 
giant clump, which are in the fast evolution phase with mass-loss, the derived masses are 
probably the upper limits. 

Stellar radii were determined using the calibration given in Alonso et al. (2000). Stellar 
ages are usually estimated with the aid of the theoretical stellar isochrones. 

In the case of the PSU-TCfA 18 star, we have measured equivalent widths of up to 195 
Fe I and 11 Fe II lines for further analysis. For this star, the basic physical parameters 
are T eff = 4246if, log .9= 2.43, v t = litems' 1 , [Fe/H] = 0.11, M = 5.5M Q , and R = 
27.4i? Q . Intrinsic uncertainties of our determinations are cT e // = 29K, erlog g = 0.09, 
avt = O.lSkms- 1 and a[Fe/H] = 0.06. 



4. V sin i measurements 

Rotation periods represent a parameter of particular importance in searches for plan- 
etary companions to red giants. As the rotation periods of these stars are very similar to 
expected orbital periods, their knowledge is critical for an unambiguous interpretation 
of observations. Any correlation of stellar activity indicators variations with the rotation 
period make a substellar companion hypothesis unlikely To estimate rotation periods 
from our spectra we have used the cross-correlation technique, as described in Benz & 
Mayor (1984). 

We have computed the CCFs by cross correlating the high S/N blue spectra with 
a numerical mask. To measure Usinz we have worked out a Usini calibration for the 
HET/HRS. To determine the do vs. (B — V) relationship, we have used 16 slow rotators 
with accurately known projected rotational velocities, preferably from Gray (1989), Fekel 
(1997) and de Medeiros & Mayor (1999). For these stars we have determined Co using the 
formula from Benz & Mayor (1984) (V sini = A^j 'cr 2 bs — a 2 ) and assuming the constant 
A = 1.9 following the Qucloz et al. (1998) and Melo et al. (2001). We have carried out a 
least-squares fit to the data of the analytical function cr = a 2 (B — V) 2 + a\ (B — V) + a , 
which yields the following calibration: cr = 15.592(5 - V) 2 - 26.753(5 - V) + 14.559. 

Using this calibration we have obtained Usini = (3 ± l)fcms _1 for our star. Adopting 
the radius for this star as determined above, we have estimated its rotational period to 
be 220 - 950 days. The large uncertainty in the rotation period is caused by uncertainties 
related to the determination of the radius and Vsini. 
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5. Stellar activity indicators 

One of possible sources of the observed RV variations in GK-giants is due to their 
pulsations. Therefore, photometric data that span long periods of time are needed for 
the interpretation of the results of our survey. Because we do not conduct our own 
parallel photometric observations, we must rely on the existing photometric databases 
like Hipparcos or NSVS (Wozniak et al. 2004). These moderate quality data provide 
time-series, which are long enough to be useful in searches for long-term periodicities. 
However, these measurements were usually performed many years prior to our RV survey. 
In the particular case of the PSU-TCFA 18 star, no detectable variability is present in 
the existing photometric data. 

5.1. Line bisectors 

The basic tool to study the origin of RV variations derived from the stellar spectra is the 
analysis of the shapes of spectral lines via line bisectors (Gray 1983). 

We have computed line bisectors for 5 strong, unblended spectral features of a moderate 
intensity, which were located close to the center of echelle orders: Cr I 663.003 nm, Ni I 
664.638 nm, Ca I 671.77 nm, Fe I 675.02 nm, and Ni I 676.784 nm. All these lines show 
well defined bisectors. 

The changes in the spectral line bisectors were quantified using the bisector velocity 
span (BVS) parameter, which is simply the velocity difference between the upper and 
the lower points of the line bisector (BVS = V3 — V2), and the bisector curvature (BC), 
which is the difference of the velocity span of the upper half of the bisector and its lower 
half (BC = (vs — V2) — (v2—vi)). It is important to examine both BVS and BC, because 
it is possible for a star to show variations in one of these parameters only. In choosing 
the span points, it is important to avoid the wings and cores of the spectral line, where 
the errors in the bisector measurements are large. For our span measurements we chose 
Vi = 0.29, v 2 = 0.57, and v 3 — 0.79 in terms of the line depth at the line core. Using 
the bisector measurements of all 5 spectral lines we have computed the average velocity 
span and curvature after subtracting the mean value for each spectral line. 

In Figure 1, we present the mean bisector velocity span (MBVS)a,nd the mean bisector 
curvature (MBC) for our star, as a function of RV. Uncertainties in the derived values of 
MBVS and MBC were estimated as standard deviations of the mean. The correlation 
coefficients were found to be r = -0.32 ± 0.10 for MBVS and r = 0.07 ± 0.02 for 
MBC. It is clear that they are not correlated with radial velocities which supports the 
planetary mass companion hypothesis. The available RV measurements are not uniformly 
distributed over the estimated period, which is visible in Fig. 1 as a scatter varying with 
RV. More observations are needed to confirm the apparent lack of correlation. 

5.2. Ha variations 

Since our spectra start at 407.6 nm we cannot investigate the variation of the Ca II K 
emission line (393.4 nm). Also the infrared Call triplet lines 849.8-854.2 are outside the 
range of our spectra. Thcrefoer, we use Ha line (656.28 nm) as a chromospheric activity 
indicator. The EW measurements of the Ha line can be made in our spectra with a 
typical precision of a few percent. In the case of our star, the measurements give a mean 
value of 1049 ± 20mA. The rms value of 20mA corresponds to 2 % variation in the EW. 
In Figure 1 we present EW measurements for Ha as a function of RV. The correlation 
coefficient of r = 0.44 ± 0.10 shows marginal relationship (probably resulting from the 
non- uniform RV coverage) which again supports the planetary hypothesis. 
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Figure 1. MBVS, MBC and EW of the Ha for PSU-TCfA 18 as a function of radial velocity. 
6. Conclusions 

A detailed knowledge of stellar parametres of red giants is very important for interpre- 
tation of their RV variations. To rule out stellar activity as the source of such variations, 
one needs precise rotation periods and several other indicators to be measured at many 
epochs. 
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